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APPROXIMATE  CALCULATION  OF  NONSTATICNARY  TEMPERATURE 
REGIME  OF  A ROCK-ICE  DAM 


Engineer  N.  V.  Ukhova 

Transactions  of  Coordinated  Meetings  on  Hydraulic 
Engineering.  Issue  No.  23,  MISI  imeni  V.  V.  Kuybyshev 


In  order  to  determine  the  possibility  of  building  a frozen  dam 
made  of  local  materials  under  certain  harsh  climatic  conditions, 
working  out  its  construction  in  detail,  and  concrolling  its 
thermal  regime,  it  is  necessary  to  be  able  to  predict  the  tem- 
perature field  of  the  dam  and  foundation.  Changes  in  the  tem- 
perature of  earth  dams  as  a function  of  the  thickness  of  perma- 
frost soil  were  investigated  by  Professor  P.  A.  Bogoslovskiy 
[1,  2,  3]. 

In  many  instances,  it  is  economically  advantageous  to  build 
a rock-fill  dam.  In  1961,  Gidroproyekt  proposed  construction 
of  rock-fill  dams  in  regions  where  permafrost  soil  occurs, 
particularly  dams  with  an  ice  core.  The  operating  regime  of 
such  dams  differs  somewhat  from  that  of  earth-fill  dams. 

During  1961-1964,  on  order  from  Gidroproyekt,  the  Department 
of  Soil  Mechanics  and  Foundations  of  MISI  imeni  V.  V.  Kuyby- 
shev, headed  by  Professor  and  Corresponding  Member  of  the 
Academy  of  Sciences  of  the  USSR  N.  A.  Tsitovich,  studied  the 
principles  of  calculation  and  planning  of  rock-ice  dams.  The 
studies  took  two  directions: 

1)  Development  of  principles  for  calculating  nonstationary 
temperature  regimes  of  rock-ice  dams; 

2)  Evaluating  the  stress-deformation  state  of  the  ice  core 
of  a rock-ice  dam. 

Some  of  the  results  of  these  investigations  have  been  published 
previously  [4,  5,  6].  The  present  article  is  devoted  to  the 
studies  which  were  conducted  in  the  first  of  the  above  direc- 
tions. The  studies  in  the  second  direction  are  set  forth  in 
another  article  in  the  present  collection. 

As  a result  of  four  years'  work,  an  approximate  method  con- 
venient tor  practical  use  was  worked  out  for  calculating  the 
nonstationary  temperature  field  of  a rock-ice  dam  and  its 
foundation.  The  results  of  the  calculations  performed  by  this 
mr  ->d  for  a concrete  dam  on  the  Tsypa  River  were  checked  on 
t,  "~'A  9/60  integrator  and  by  a series  of  experiments  on  a 
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physical  model.  Comparison  of  the  calculated  and  experimental 
temperature  fields  showed  good  agreement.  This  indicated  that 
this  calculation  method  was  suitable  for  practical  use. 

To  develop  a method  of  calculating  the  nonstationary  tempera- 
ture regime  of  a rock-ice  dam,  the  following  diagram  was  used 
(Figure  1) . 


Figure  1.  Design  Diagram  for  Rock 

Dam.  I,  Rock  fill  with  ice;  II,  Ice  core; 

III,  Dry  rock  fill;  IV,  Foundation. 


The  nonstationary  temperature  regime  of  a rock-ice  dam  and 
its  foundation  must  be  viewed  in  the  general  case  as  a spa- 
tial problem.  At  the  present  time,  this  problem  cannot  be 
solved  in  practice.  If  we  exclude  from  consideration  the 
range  in  which  there  is  contact  between  the  dam  and  the 
edges  of  the  canyon,  the  problem  may  be  considered  a plane 
one.  However,  even  this  problem  cannot  be  solved  in  a closed 
form  due  to  its  complexity.  We  must  deal  with  a temperature 
field  which  changes  as  a function  of  time  wxthin  the  specific, 
limited,  and  finite  dimensions  of  an  area  with  different 
thermal  physical  characteristics  and  different  initial  and 
boundary  conditions  in  separate  sections.  Therefore,  when 
working  out  the  method  of  calculation,  the  following  assump- 
tions were  made. 

1.  The  general  plane  problem  is  expanded  into  a number  of 
linear  problems  and  one  plane  problem,  each  of  which  is  solved 
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independently  of  the  others.  The  general  temperature  field 
ct  any  moment  in  time  is  composed  of  partial  temperature 
fields.  The  linear  problems  are  solved  on  t.he  basis  of 
classical  problems  of  the  theory  of  thermal  conductivity 
tor  semi-limited  rods  with  lateral  insulation  [7J.  The 
plane  problem  is  the  distribution  of  temperature  in  the  foun- 
dation, and  is  solved  by  the  method  of  finite  differences  11]. 

2.  The  temperature  distribution  in  the  general  case  is  stud- 
ied by  breaking  up  the  nonstationary  process  into  a series 

of  stationary  processes,  alternating  with  one  another. 

3.  The  construction  process  is  not  taken  into  account.  The 
time  is  counted  from  when  the  dam  has  already  been  built  and 
the  reservoir  is  filled.  At  this  moment,  certain  boundary  con- 
ditions are  established  which  remain  constant  during  the  ex- 
istence of  the  structure. 

4.  The  spaces  in  the  fill  at  the  upper  wedge  are  filled  with 
ice.  The  initial  temperature  of  all  the  zones  in  the  body  of 
the  dam  is  0®. 

5.  In  the  base  of  the  dam,  there  is  a layer  with  no  cracks, 
whose  cooling  and  heating  take  place  without  any  phase  shifts. 

Then  the  temperature  regime  of  the  rock-ice  dam  and  its  foun- 
dation (Figure  1)  is  solved  on  the  basis  of  the  following 
problems : 

1)  Thawing  of  the  upper  wedge  of  the  dam  from  the  heat  in 
the  reservoir.  Section  I-I; 

2)  Cooling  of  parts  of  the  upper  and  lower  wedges  and  the 
crest  of  the  dam  due  to  the  negative  air  temperature.  Section 

3)  Heating  of  the  foundation  at  a considerable  distance  from 
the  upper  slope  by  the  heat  from  the  reservoir,  Section  III- 
III; 

4)  Cooling  of  the  foundation  at  a considerable  distance  from 
the  lower  slope  due  to  the  negative  temperature  of  the  air. 
Section  IV-IV; 

5)  Thawing,  heating  or  cooling  of  the  cushion  of  the  dam  under 
the  influence  of  the  temperature  field  of  the  foundation.  Sec- 
tion V-V; 
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6)  Thawing,  heating,  or  cooling  of  the  laminar  layers  (body 
of  the  dam  - core  - foundation).  Section  VI-VI; 


7)  Taking  into  account  sets  of  transformations  at  the  boun- 
dary between  the  frozen  and  thawed  zones; 

8)  Calculation  of  forced  cooling  (cooling  galleries); 

9)  Calculation  of  the  temperature  field  of  the  foundation. 

The  boundary  conditions  are  general  for  all  the  problems 
(Figure  1) : 


t,  (t)  = const  >0 
f:  (■*)  “■  const  <"  0 
.•,(t«0)  -r.(T-o)-o 
/,(■:)  Ai  const  < 0 


(1) 


where  t = f(x)  is  the  temperature  of  the  foundation,  x being 
the  depth  from  the  cushion  of  the  dam.  In  these  and  subse- 
quent formulas,  t is  the  calculation  time  interval  and  t is 
the  temperature. 


In  the  following  we  shall  examine  the  solutions  of  the  above 
problems. 


Problem  1.  Thawing  of  the  upper  wedge  due  to  the  heat  from 
the  reservoir  is  accompanied  by  thawing  of  the  ice  in  the 
spaces  in  the  rock  fill.  Since  = 0®,  when  the  slope  is 

rather  long,  the  depth  of  the  thawing  at  a moment  in  time  t, 
calculated  along  a perpendicular  to  the  slope,  will  be  deter- 
mined by  the  Stefan  formula: 


(2) 


where  is  the  coefficient  of  thermal  conductivity  of  the 
thawed  zone; 

p is  the  latent  heat  of  melting  of  the  ice; 

Y is  the  specific  gravity  of  the  skeleton  of  the  rock 
fill; 

w is  the  moisture  content  of  the  rock  fill  in  fractions 
of  a unit. 


With  a certain  granulometric  composition  for  the  fill,  a con- 
siderable influence  on  the  rate  of  thawing  of  the  upper  wedge 
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may  be  exerted  by  convective  heat  exchange  processes.  A 
method  was  worked  out  to  calculate  the  influence  of  convec- 
tive heat  transfer  within  general  heat  exchange,  with  thaw- 
ing of  the  upper  wedge.  This  took  two  forms:  the  first  in- 
volved joint  solution  of  heat  transfer  equations  by  convec- 
tion and  thermal  conductivity  as  they  apply  to  a concrete 
scheme,  similar  to  the  one  used  by  other  researchers  [81. 

The  second  involved  determining  the  speed  (or  flow) , and 
from  this  the  heat  flux,  with  filtration  of  water  in  the 
granular  material  under  the  influence  of  the  difference  in 
hydrostatic  pressure  caused  by  the  different  density  of  the 
water  at  different  temperatures.  Since  the  heat  gradient  is 
very  small  with  natural  convection,  the  determination  of 
analytical  relationships  is  based  on  Darcy's  Law. 

Both  of  these  approaches,  based  upon  different  original  con- 
siderations, lead  to  the  same  result.  Unfortunately,  due  to 
the  limited  scope  of  this  article,  we  cannot  outline  in 
detail  the  method  of  calculation  of  convection  and  we  are 
forced  to  present  only  the  final  result:! 


where 


•1  - 1 = O.OtHM) 


htj  - t,)  n-^n'Cph  , 


(3) 


j—  is  the  coefficient  of  the  increase  in  thermal 
T conductivity  due  to  the  presence  of  convective 
currents  (X,,  is  the  coefficient  of  the  effec- 
tive thermal  conductivity,  assuming  that  all  the 
heat  transmitted  by  contact,  in  other  words  by 
convection  and  thermal  conductivity,  is  trans- 
mitted only  by  thermal  conductivity) ; 


0.0006  is  the  coefficient  which  takes  into  account  the 
hydraulic  resistance  of  the  granular  layer  and 
the  nature  of  the  convective  currents  in  the  fill 
[8]; 


S is  the  coefficient  of  volumetric  expansion  of 
water; 

ta  is  the  temperature  of  water  at  maximum  density; 
ti  is  the  temperature  of  water  when  frozen; 
g is  the  acceleration  due  to  gravity; 

Y is  the  density  of  water; 
n is  the  void  coefficient  of  the  fill; 


This  question  is  discussed  in  detail  in 


[6]. 
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^3 
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is  the  specific  thermal  capacity  of  the  water; 

is  the  dynamic  coefficient  of  viscosity; 
is  the  coefficient  of  thermal  conductivity  of  the  fill; 
is  the  equivalent  diameter  of  the  rock  fill,  deter- 
mined by  the  formula  of  Ye.  A.  Zamarin; 
is  the  depth  from  the  waterline  in  the  reservoir  at 
which  the  influence  of  convection  is  taken  into  ac- 
count. 


When  establishing  the  constant  parameters  characterizing  the 
water  and  the  material  of  the  medium  (B,  tj,  ti,  y,  g,  Cp»  y, 
X)  Formula  (3)  becomes 


?-l  - {3') 

The  solution  to  (3')  is  in  the  form  of  a graph  (Figure  2). 

On  this  graph,  after  determining  the  value  of  d^  for  the 
fill,  with  known  values  for  n and  h,  we  can  determine  - 1. 
There  are  three  areas  on  the  graph:  (1)  when  O - 1 < 0.1 
convection  has  practically  no  effect;  (2)  when  0.1<'9-1<3 
convection  affects  the  positions  of  the  isotherms  and  must  be 
tciken  into  account;  (3)  when  3 < <J)  - 1 convection  predominates 
over  thermal  conductivity. 

Depending  on  the  values  obtained  for  (J)  - 1,  convection  must 
be  taken  into  account.  When  it  is  necessary  to  take  convec- 
tion into  account,  after  determining  't>  - 1 and  then  using 
this  value  to  determine  X3  = 4>A<p,  we  can  substitute  X in  (2) 
instead  of  X.p. 

Problem  2.  The  cooling  of  portions  of  the  upper  and  lower 
wedges  and  crest  of  the  dam  is  solved  on  the  basis  of  the  prob- 
lem of  cooling  of  semi-bounded  rods  with  lateral  insulation 
[7]. 

The  solution  has  the  following  form: 


^(J:.  = 1 (4) 

V j 1 (It  } 


where  (x,  t)  is  the  temperature  at  a point  at  a distance  x 
along  the  normal  to  the  surface,  from  which 
cooling  is  determined  at  a moment  t; 
a is  the  coefficient  of  thermal  conductivity  of 
the  material  in  the  zone  in  question  in  a frozen 
state . 
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Figure  2.  Graph  for  Determining  Convection  in 
General  Heat  Exchange  as  a Function  of  Porosity 
and  Effective  Rock  Diameter.  On  the  diagram: 

(A)  Region  I - convection  has  practically  no 
effect  and  may  be  disregarded.  (B)  Region  II  - 
convection  has  a considerable  influence  on  the 
position  of  the  isotherms  and  must  be  taken  into 
account.  (C)  Region  III  - convection  predominates, 
and  thermal  conductivity  may  be  disregarded. 
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Problem  3.  Heating  of  the  foundation  at  a considerable  dis- 
tance  from  the  upper  slope  on  the  basis  of  the  problem  of 
Academician  A.  V.  Lykov  [7]  on  the  heating  of  a semi-bounded 
rod  with  lateral  insulation.  The  general  solution  has  the 
form: 


{ I / ,;)  ftp 

I 

/ ('•  • M./. 


4ut 


(5) 


where  t{x,  t)  is  the  initial  temperature  at  a point  with  co- 
ordinate X at  a moment  in  time  t ; 
a is  the  coefficient  of  thermal  conductivity  of 
the  material  of  the  foundation; 
erfc  n = 1 - erf  n,  where  erf  n is  the  integral  of  the 
Gaussian  errors; 

fiO  is  a function  which  determines  the  initial  dis- 
tribution of  the  temperature  in  the  foundation. 

When  f(C)  = ti  = const  (1)  becomes; 


/i.v,  -) 


[t,-  ;,)eric----r=- 
2 I ir 


(5‘) 


Problem  4.  Cooling  of  the  foundation  at  a considerable  dis- 
tance  from  the  lower  slope  xs  solved  analo^usly,  by  svobstituting  t; 
for  tz  in  (5).  Otherwise,  the  method  of  solving  this  prob- 
lem is  the  same. 

Problem  5.  Thawing,  heating,  or  cooling  of  the  cushion  of 
the  dam  under  the  influence  of  the  temperature  field  in  the 
foundation  must  be  viewed  from  two  angles; 

1)  The  foundation  of  the  dam  is  a talik,  in  other  words  the 
rock  is  at  a positive  temperature; 

2)  The  base  of  the  dam  is  made  of  permafrost  soil  with  a 
finite  thickness  (hm) , with  the  rock  beneath  at  a positive 
temperature.  The  solutions  to  these  problems  are  found  on 
the  basis  of  the  heat  balance  equations  proposed  by  I.  S. 
Moiseyev  [5,  6] . In  the  first  case  the  thawing  of  the  fro- 
zen upper  wedge  on  the  base  side  will  be  expressed  by  the 
formula: 


8 


The  heating  of  the  dry  lower  wedge  on  the  base  side  is  des- 
cribed by  the  formula: 


I 


•*•1  = -f  .• 


(8) 


In  these  formulas,  X with  the  corresponding  subscripts  repre- 
sents the  coefficients  of  thermal  conductivity  of  the  material 
in  the  zones  of  the  dam  [in  (6) , taking  convection  into  ac- 
count if  necessary) ; and  Cq  are  the  specific  heat  capacity 
of  the  dry  fill  and  water;  tj^  is  the  temperature  at  the  con- 
tact between  the  base  and  dam;  hi  is  the  depth  at  which  the 
temperature  remains  constant  during  the  calculation  period 
(according  to  our  calculations  on  the  EGDA  9/60,  it  may  be 
assumed  to  be  300-500  meters) ; p is  the  specific  heat  of  melt- 
ing of  ice  and  w is  the  moisture  content  by  weight  of  the 
fill. 


in  the  second  case  the  cooling  of  the  frozen  upper  wedge  on 
the  base  side  is  expressed  by  the  formula 


).,(-/.)  At 


(6') 


The  cooling  of  the  ice  core  from  the  base  side  is  expressed 
by; 


'*‘.1  ( ^«) 
'•  2 
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(7') 


w 


The  cooling  of  the  dry  lower  wedge  from  the  base  side  is 
described  by: 


*•«(-/.)  At 


(8*) 


The  movement  of  the  lower  limit  of  the  permafrost  in  the 
foundation  C:  during  a brief  period  of  time  (3-5  years)  from 
the  moment  the  dam  goes  into  service  is  expressed  by  the 
formula: 


The  symbols  in  these  formulas  are  the  same  as  those  used 
before,  but  the  subscript  "M"  shows  that  the  characteristics 
were  plotted  for  material  in  a frozen  state. 

Problem  6 . The  thawing,  heating,  or  cooling  of  the  laminated 
thickness  is  considered  in  those  cases  when  the  zero  isotherm 
runs  from  an  area  with  homogeneous  thermal  physical  charac- 
teristics into  an  area  with  different  ones.  It  is  necessary 
in  this  connection  to  use  the  recommendations  of  V.  S.  Lukyanov 
and  M.  D.  Golovko  regarding  the  average  derived  characteristics 
[9]  . 

Problem  7.  Calculation  of  the  aggregate  conversions  takes 
place  when  in  some  region  (the  core  for  example)  there  is  a 
frozen  zone  and  a thawed  zone,  with  an  unstable  temperature 
regime  and  constant  temperatures  at  the  boundaries  of  these 
zones.  Then  the  shifting  of  the  zero  isotherm  and  the  mixing 
of  the  remaining  isotherms  associated  with  it  may  be  considered 
on  the  basis  of  the  known  condition  at  the  interface  [7] : 

A;  - 1 < « J^rdd  | | yrad  t,  \ — . (10) 

where  ) grad  tj^j)  and  | grad  trji|  are  the  temperature  gradients 
in  the  frozen  and  thawed  zones. 

The  positive  value  of  AC  indicates  an  increase  in  the  frozen 
zone.  Formula  (10)  was  compiled  for  an  ice  core.  For  frozen 


fill,  the  term  on  the  right-hand  side  of  this  formula  which 
takes  into  account  the  aggregate  conversions  will  be  as 
follows: 


At 


Problem  8.  Calculation  of  forced  cooling  is  considered  for 
a horizontal  gallery  with  a diameter  of  2ro  with  a constantly 
operating  source  of  cold  (runoff)  with  a thickness  of  i-  qg) • 

The  solution  ’;as  obtained  on  the  basis  of  the  problem  of 
Academician  A.  V.  Lykov  on  temperature  distribution  in  two 
contiguous  semi-limited  rods  with  lateral  insulation  from  a 
constantly  moving  source  of  cold  at  the  boundary  between  them 
[7].  The  solution  has  the  following  form: 


M.c,  r)  = -• 


2</..  \ ar‘,~  ^«-i) 


/eric 


X 


2 ) /J.T 


(11) 


The  symbols  in  the  formula  are  the  same  as  those  in  the  formu- 
las given  above. 


Use  of  (11)  to  plot  the  temperature  field  in  a system  as  com- 
plicated as  an  inhomogeneous  laminated  medium  (different  zones 
of  the  dam  and  foundation)  is  only  possible  with  the  following 
assumptions : 

1.  Temperature  is  determined  as  a function  of  time  along  a 
radial  ray  extending  from  the  center  of  the  gallery,  and  not 
from  its  surface;  in  other  words,  it  is  considered  that  the 
source  of  the  cold  is  a point  source  and  the  entire  gallery 
is  filled  with  ice.  This  reduces  the  effect  of  the  influence 
of  the  gallery,  producing  a reserve  of  thermal  strength. 

2.  In  accordance  with  the  second  assumption  set  forth  at  the 
beginning  of  the  article,  the  calculation  is  performed  for 
certain  moments  in  time  Ti,  t 2 , T3,  ...  from  the  beginning 

of  reckoning,  assuming  a steady  state  condition  between  these 
moments.  Then  for  each  moment  in  time,  instead  of 
we  substitute  in  (11)  Ti  - 0;  ta  - t 1 ; T3  - Ta;  etc.  The 
initial  temperature  tg  at  each  moment  in  time  is  selected  to 
be  an  average  within  the  radial  ray  along  which  temperature 
distribution  is  calculated. 
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The  thickness  of  th<;  cold  source  (-  qo)  is  expressed  as  a 
function  of  the  method  of  cooling  used  and  the  cooling  equip- 
ment. Approximate  calculation  of  this  factor  must  be  based 
upon  an  analysis  of  the  thermal  balance  of  the  thawing  zones. 

Problem  9.  In  solving  the  temperature  field  in  the  foundation, 
by  contrast  with  the  previous  one-dimensional  problems,  a 
two-dimensional  one  is  considered.  For  this  purpose,  an 
approach  is  used  which  was  developed  by  P.  A.  Bogoslovskiy 
and  is  presented  in  his  paper  [1] . 

The  temperature  regimes  of  concrete  structures  can  be  calcu- 
lated conveniently  under  certain  conditions  on  the  basis  of 
the  problems  set  forth  above,  as  follows: 

1.  Using  the  structural  diagram  of  the  structure,  the  phys- 
ical, mechanical,  and  thermal  physical  characteristics  of 
various  structural  zones  of  the  body  of  the  dam  and  founda- 
tion are  determined.  In  individual  zones  it  is  necessary  to 
determine  only  those  regions  where  the  characteristics  differ 
considerably  from  one  another:  n by  20%;  y by  0.2  m/m^; 

a.  A,  C by  30  to  50%.  In  the  remaining  cases,  the  average  de- 
rived characteristics  of  the  zones  are  established.  On  the 
basis  of  comparison  and  analysis  of  these  characteristics, 
the  design  diagram  of  the  dam  and  foundation  are  then  prepared. 

2.  The  boundary  conditions  are  established  for  the  design 
diagram. 

a)  Geometric  dimensions  corresponding  to  the  dimensions  on 
the  design  diagram; 

b)  The  initial  conditions  = to  = 0 for  temperature  dis- 
tribution in  the  foundation  f(C); 

c)  Boundary  conditions  ti,  ta. 

3.  The  influence  of  convection  on  the  thawing  of  the  upper 
wedge  is  estimated: 

a)  By  using  the  formula  of  Ye.  A.  Zamarin,  d^  is  determined; 

b)  Using  the  graph  in  Figure  2,  ({>  - 1 and  A^  are  determined. 

4.  The  temperature  field  of  the  dam  and  foundation  is  built 
up  at  certain  moments  in  time  from  the  arbitrary  beginning  of 
calculation.  For  average  and  large  dams,  it  is  recommended 
that  the  following  moments  in  time  be  selected:  Xo  = 0 years. 
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Ti  = 5 years;  T2  = 10  years,  T3  =20  years,  T4  = 40  years; 

Ts  = 75  years,  Te  = 100  years. 

5.  In  the  design  diagram  for  each  t^,  cross  sections  are 
selected  along  which  temperature  calculations  may  be  performed. 

The  cross  sections  must  reflect  as  fully  as  possible  the  tem- 
perature changes  in  the  areas  of  the  dam  and  foundation  in 
question. 

b.  The  temperature  values  are  calculated  for  , using  the 
selected  cross  sections  without  a cooling  gallery.  On  the 
basis  of  these  values,  a general  temperature  field  is  con- 
structed. 

7.  For  T2,  cross  sections  are  again  selected  and  the  tempera- 
ture and  structure  of  the  general  temperature  field  are  then 
calculated  once  again. 

8.  As  soon  as  the  temperature  fields  have  been  calculated 
for  the  entire  design  time  period,  the  thermal  stability  of 
the  dam  and  foundation  are  analyzed  and  a decision  is  made 
regarding  the  need  for  forced  cooling.  If  forced  cooling  is 
necessary,  the  influence  of  the  gallery  is  superimposed  upon 
the  temperature  field  calculated  at  xj,  T2,  x^,  etc.  by  (11). 

Figure  3 shows  an  example  of  the  temperature  field  of  a dam 
and  foundation  constructed  according  to  the  proposed  method 
for  the  type  of  rock-ice  dam  planned  by  Gidroproyekt  for  the 
Tsypa  River,  with  x = 10  and  75  years. 

I 
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